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Abstract-Vlz  have  developed  a  generic  ionic  model  of  single-cell 
rabbit  sinoatrial  node  (SAN)  electrical  activity,  employing 
Markov-state  kinetics  for  membrane  ionic  currents.  By  suitable 
choice  of  parameters,  the  model  is  able  to  accurately  reproduce 
a  wide  range  of  action  potential  waveforms  recorded  from 
individual  SAN  cells  in  the  laboratory.  Advantages  of  this 
approach  include  the  ability  to  quantitatively  investigate  the 
contribution  of  ion  currents  underlying  pacemaker  and  action 
potential  activity  from  any  given  SAN  myocyte  recording. 
Keywords  -  Sinoatrial  node,  action  potential,  mathematical 
model. 


I.  Introduction 

Mammalian  hearts  generate  self-sustaining  rhythmic 
activity  from  trigger  impulses  arising  from  a  region  of 
pacemaker  cells  located  in  the  wall  of  the  right  atrium  known 
as  the  sinoatrial  node  (SAN).  Biophysically-detailed 
mathematical  models  of  single  cell  SAN  electrical  activity 
have  been  extensively  published  [1,  2],  all  based  on  Hodgkin- 
Huxley-type  formalism  to  simulate  underlying  ionic 
membrane  currents.  Although  yielding  generic  action 
potential  waveforms  crudely  similar  to  those  recorded  in  the 
laboratory,  none  of  the  models  were  able  to  accurately 
reproduce  action  potential  waveforms  from  any  specific  SAN 
cell.  This  limits  the  usefulness  of  such  models  to 
quantitatively  assess  the  contribution  of  underlying 
membrane  currents  to  pacemaker  activity,  as  well  as  the 
effects  of  various  pharmacological  agents,  in  any  given  SAN 
myocyte. 

There  is  a  need  therefore,  for  a  cell-specific  approach  in 
modelling  SAN  pacemaker  activity,  which  can  accurately 
simulate  action  potential  waveforms  recorded  in  the 
laboratory.  In  this  study,  we  describe  a  new  single  cell  SAN 
model  employing  Markov-state  kinetic  descriptions  of  ion 
currents,  capable  of  highly-accurate  reproduction  of  action 
potential  waveforms  recorded  from  specific  SAN  myocytes. 


II.  Model  Description 

The  SAN  cell  model  developed  in  this  study  consists  of  12 
membrane  currents,  acting  in  parallel  across  a  capacitive  cell 
membrane,  with  transmembrane  E  potential  given  by 
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—it. 


dt  C 


(i) 


with  total  membrane  current  itot  given  by  the  algebraic  sum  of 
its  12  constituents: 


hot  -  lf  +  i K,r  +  lK,s  +  h,Na  +  *  Na  +  %  +  - 

U) 

lNaCa  +  ho  +lCa,L  +lCa,T  +lb,K  +lb,Cl 

A  brief  description  of  the  currents  in  (2)  is  summarized  in 
table  1 .  All  time-dependent  currents  were  described  using  3 
or  4  Markov-state  kinetic  schemes,  as  outlined  in  table  1 . 
These  kinetic  schemes  represent  conformational  changes  in 
membrane-spanning  proteins  controlling  the  flow  of  ions 
across  the  cell  membrane.  With  reference  to  table  1,  stable 
states  of  the  underlying  channel  are  denoted  by  C  (closed),  A 
(activated  or  open),  and  Ib  I2  (inactivated).  Forward  and 
reverse  rate  constants  between  states  is  a  function  of 
membrane  potential  according  to 


k  =  k^+- 


k0  -kx 


1  +  exp 


(E-E  5p) 

E  slope 


(3). 


This  equation  defines  a  sigmoidal-type  profile  for  rate  A;  as  a 
function  of  membrane  potential  defined  by  4  parameters:  k ^ 
k0  (the  values  of  k  at  £=±00),  E50  (the  E  value  at  50%  of  Us 
range)  and  Eshpe  (related  to  the  slope  at  E  =  E5q).  Each  rate  in 
the  kinetic  schemes  of  table  1  is  given  by  (3),  each  having  its 
own  set  of  4  parameters. 

The  general  form  for  all  time-dependent  currents  was 
given  by 

i  =  gA(E-Erev)  (4) 


with  g  denoting  the  membrane  conductance,  Erev  the  reversal 
potential,  and  A  the  proportion  of  membrane  channels  in  the 
open  state.  Default  conductances  and  rate  parameters  in  (3) 
and  (4)  were  chosen  to  fit  published  voltage-clamp  data  for 
each  ion  current. 

Using  the  default  values  of  all  model  parameters  as  a 
starting  point,  a  total  set  of  175  parameters  were  adjusted  to 
fit  a  variety  of  pre-recorded  single  cell  action  potential  SAN 
waveforms  using  a  custom  non-linear  optimization  routine 
developed  for  this  purpose.  All  rate  parameters  were 
constrained  to  lie  within  ±30%  of  default  values,  and  all  ion 
current  peaks  were  constrained  to  lie  within  physiological 
limits.  This  ensured  that  the  contribution  of  all  membrane 
currents  to  pacemaker  activity  was  in  accordance  with 
commonly  reported  experimental  values.  A  final  constraint 
was  that  all  ionic  species  K+,  Na+,  Ca2+  and  Cl"  should  exhibit 
no  cycle-cycle  drift. 
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III.  Results 


Table i 


SUMMARY  OF  MEMBRANE  CURRENTS  IN  SAN  CELL  MODEL 


Symbol 

Name 

Kinetic  Scheme 

Hyperpolarization-activated  current 

A 

c^±  c 

2  “3  1 

lK,r 

Rapid  delayed-rectifying  K+  current 

a, 

A  - **  1, 

tii.  h 

f  '  a4  ^2 

lK,s 

Slow  delayed-rectifying  K+  current 

A 

c  ^=±  c 

2  «3  1 

h.Na 

Background  Na+  current 

(time-independent) 

iNa 

TTX  sensitive  Na+  current 

a, 

A  — **  1, 

fli.  I' 

h  '  a4  ^2 

ip 

Na+-K+  pump  current 

(time-independent) 

iNaCa 

Na+-Ca2+  exchange  current 

(time-independent) 

ho 

Transient  outward  current 
(incorporating  4-AP-sensitive 
sustained  outward  current  isus) 

a, 

A  — **  1, 

fli.  1- 

h  '  a4  ^2 

ica,L 

L-type  Ca2+  current 

A  , 

•K+ 

i2 

ica.T 

Transient  Ca2+  current 

A 

c^±  c 

2  «3  1 

h,K 

Background  K+  current 

(time-independent) 

h.ci 

Background  Cl'  current 

(time-independent) 

Fig.  1 .  Computational  fit  obtained  to  isolated  single  cell  sinoatrial  node 
recording.  Top  trace  illustrates  the  near-exact  overlap  in  action  potential 
between  the  real  cell  and  an  ionic  model.  Lower  traces  reveal  underlying  model 
membrane  currents.  All  initial  variables,  parameters  and  peak  currents  of  the 
model  were  constrained  to  lie  within  physiological  limits.  Membrane  currents 
were  also  constrained  to  yield  zero  mean  ionic  fluxes  per  cycle  for  steady-state 
periodic  activity. 

Figure  1  illustrates  an  example  of  one  model  fit  to  action 
potential  waveform  data  recorded  from  an  enzymatically- 
isolated  SAN  cell  using  the  Nystatin  patch  technique.  RMS 
error  between  the  model  waveform  and  the  smoothed 
experimental  recording  was  approximately  1  mV.  The  results 
of  figure  1  reveal  a  marked  improvement  over  previous 
models  in  matching  action  potential  data  recorded  from  real 
SAN  cells.  In  addition,  underlying  current  waveforms  and 
amplitudes  shown  in  this  figure  are  well  in  agreement  with 
generally  accepted  experimental  findings. 

IV.  Conclusion 

Ionic  mechanisms  underlying  pacemaker  activity  in  SAN 
cells  can  be  investigated  using  the  computational  models  of 
cellular  electrical  activity  developed  in  this  study.  A  major 
improvement  over  existing  models  is  that  model  parameters 
have  been  adjusted  to  accurately  reproduce  action  potential 
waveforms  recorded  from  specific  cells.  This  study  represents 
the  first  stage  in  highly  accurate  modelling  of  ionic  activity  in 
any  individual  pacemaker  cell. 
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